Dependence of Ridge Formation on Trigger Azimuth: Correlated Emission Model 
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Ridge formation in near-side correlation in heavy-ion collisions is studied in the framework of a 
phenomenological model, called Correlated Emission Model (CEM). Successive soft emissions due 
to jet-medium interaction lead to the enhancement of thermal partons which follow the local flow 
directions. The correlation between the flow direction and the semihard parton direction is the 
major factor that causes the ridge formation to depend on the trigger direction relative to the 
reaction plane. With the use of a few parameters we have been able to reproduce the data on the 
ridge yields as functions of the trigger azimuthal angle for different centralities. An inside-outside 
asymmetry function is proposed to further probe the characteristics of the azimuthal correlation 
function. Insights are provided for the understanding of some detailed aspects of the centrality 
dependence. 
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I. INTRODUCTION 

Experimental studies of dihadron correlations in 
heavy-ion collisions have revealed important information 
about the interaction of jets with dense medium [1, _2, 
The observation of ridges on the near-side of triggers Q 
has further stimulated a large number of theoretical in- 
vestigations on the dynamical origin of ridge formation 
More detailed properties of the ridges are needed 
to discriminate the different models. Recently, Feng has 
reported on the dependence of the ridge yield on the az- 
imuthal angle (jjs between the trigger angle (f>T and the re- 
action plane '^rp At mid-central (20 - 60%) Au-Au 
collisions the ridge yield is found to decrease rapidly with 
increasing 0s , a feature that has not been hinted by pre- 
vious experiments, nor predicted theoretically. The aim 
of this paper is to reproduce that feature in the frame- 
work of a new model which we call Correlated Emission 
Model (CEM). 

Various properties of the ridges are already known ex- 
perimentally. The ridge yield increases significantly with 
A'part, but decreases only slightly with trigger momen- 
tum [1] . The dependence on the associated particle 
momentum p^^°'^ is exponential, its slope being nearly in- 
dependent of p^'^ The baryon/meson (B/M) ratio in 
the ridge is large, comparable to that in the inclusive dis- 
tributions [3, [3] ■ Since the above properties are mostly 
revealed at intermediate pt (< 6 GeV/c) where the anal- 
yses have been carried out, those ridges are strongly cor- 
related to semihard jets [16j] . Because of energy loss of the 
semihard parton in traversing the dense medium, most of 
those jets are created near the surface. The lost energy 
enhances the thermal partons, which give rise to the ridge 
that is exponential in p^^""^. The large B/M ratio in the 
ridge suggests that the hadronization process is recombi- 
nation [17, 18;, 19]. When averaged over all trigger angle 
(j>T, the A(j) distribution of the ridge is confined to the 
range |A0| < 1. All these properties can be incorporated 



in a suitable model. The challenge now is to reproduce 
the </)s dependence of the ridge yield as found in [l^ . 

It should be mentioned that there is another piece of 
data on the length of the ridge in Arj that is new and 
striking. PHOBOS has found that it extends over the 
range —4 < A77 < 2 [20[. The model that we describe 
in this paper aimed at treating the (f>s dependence in the 
transverse plane is not suitable for describing the A77 de- 
pendence in the longitudinal direction for the same rea- 
son that correlations at large px and large p^ involves 
different physics. We envision that the long-range Arj 
correlation is due to the interaction between high-Zcy hard 
partons and the large-fci, bulk partons, a subject that is 
not our concern in this paper. Henceforth, we shall con- 
sider only correlation in the transverse plane with jvyj < 1, 
as it is restricted in the analysis in |13| . 

The range of px studied for the (j>s dependence is lim- 
ited to 3 < < 4 GeV/c and 1.5 < < 2.0 
GeV/c [l^. In fact, the experimental evidence for the 
exponential behavior of the ridge particles does not ex- 
tend to p^^""^ < 2 GeV/c, there being no data points on 
the distribution below that threshold Q • Thus we 
shall not investigate the pt dependence in our model; to 
do so would involve issues, such as the recombination of 
thermal partons, that are not critical to our focus on the 
azimuthal problem. 

We shall assume that the trigger angle 4>t is the angle 
of the semihard parton that initiates the trigger particle. 
Event-by-event the two angles may not be identical, but 
on average that is not a bad approximation. It is im- 
portant to recognize that there is a significant difference 
in time scales between the dynamics of semihard partons 
and the formation of ridge particles. The latter takes 
place at late time and is influenced by collective flow. 
Compared to that time scale semihard scattering occurs 
almost instantaneously and is sensitive to the initial con- 
figuration of the collision system. The difference in the 
time scales also leads to a difference in the azimuthal an- 
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gles of the semihard parton (f>s and of a ridge particle (j). 
Our concern will mainly be in the correlation between (pg 
and (p in the transverse plane, with parton momentum 
kx and hadron momentum pT essentially fixed in narrow 
ranges appropriate for the experimental p^'^ and p^^°'^. 

Since the effects of semihard scattering cannot be cal- 
culated reliably from first principles, we build a model 
that incorporates all the features mentioned above con- 
cerning ridge formation. The main characteristic about 
the model is the relationship between the local flow of 
the medium and the semihard parton that traverses that 
medium. The interaction between the two leads to a cor- 
relation that increases the effect of soft emission when the 
enhanced thermal partons are all within a cone around 
the jet direction. If the energy loss to the medium due 
to successive soft emission is restricted to the vicinity of 
the semihard parton in general agreement with the direc- 
tion of collective flow, the enhancement of thermal par- 
tons that leads to ridge formation is clearly more effective 
than if the flow direction is normal to the jet direction. 
We shall show that this correlation in the jet and flow di- 
rections will not only reproduce the (ps dependence of the 
ridge yield, but also lead to other consequences that can 
readily be tested by analyzing available data in appropri- 
ate ways. One such measure is the inside-outside asym- 
metry function whose properties we can predict. Another 
result is the surprising feature that at fixed small (ps the 
ridge yield per trigger has a bump as a function of impact 
parameter, i.e., it does not decrease monotonically with 
increasing peripherality . 

The focus of our study in this paper is on the correla- 
tion between trigger and ridge particles on the near side. 
The structure of the away-side distribution involves ad- 
ditional physics not relevant on the near side because of 
trigger bias that results in higher jet momentum on the 
away side in order to penetrate the bulk medium, lead- 
ing to such phenomenon as the formation of double peak. 
Nevertheless, the physics of ridge formation considered 
here may well be pertinent to the properties of the away- 
side peaks, a subject worthy of dedicated investigation in 
a separate study. 



II. GEOMETRICAL CONSIDERATIONS OF 
SEMIHARD PARTONS 

Since the trigger momentum is in the range 3 < p^'^ < 
4 GeV/c, the semihard partons are at least also in that 
range. They lose some momentum in traversing medium, 
but gain some in recombination with thermal partons to 
form hadrons. With parton kx > i GeV/c, the time 
involved in its creation is roughly k^^ < 0.1 fm/c, which 
is short enough to be sensitive to the initial configuration 
of the system. Let the point of creation in the almond- 
shaped overlap region be labeled by P with coordinates 
(xo,yo) in transverse plane, the x-axis being in the 
reaction plane, i.e. rp = 0, and the y-axis being along 
the long side of the almond. The semihard parton can be 



scattered into any azimuthal angle, which may differ from 
the trigger angle (f)T in any event, but on average they 
may be taken to be the same. Thus hereafter we shall 
assume that the semihard parton's angle is (j)s, the same 
as measured in the experiment jlSj], i.e. (ps — (pT — ^ rp- 

Because of energy loss to the medium a semihard par- 
ton that can trigger an event is not likely to originate 
from the deep interior of the system. Neither can it be 
right at the surface if it is to generate a ridge, since it 
has to interact with the medium and lose enough energy 
to enhance the thermal partons. If P is a point at, say, 
1 fm away from the surface, then in a time interval of 
1 fm/c for the parton to reach the surface, the medium 
will have expanded and the geometry of the system may 
be better described by an expanding ellipse. There is no 
reliable way to treat the problem, since hydrodynamics 
is not applicable without strong assumptions about fast 
thermalization at time r < 1 fm/c, and QCD is not per- 
turbative for semihard scattering. In the model that we 
construct to describe the process, we consider a straight- 
line trajectory for the parton from P toward the bound- 
ary at angle (ps, ignoring the recoil parton that moves 
toward the interior and is absorbed by the bulk medium. 
The boundary is to be described by an ellipse. The dis- 
tance between P and the boundary measured along the 
trajectory is to be denoted by t (not time). At points 
along the trajectory soft emission occurs that leads to 
the development of the ridge particles, the distribution 
of which will be discussed in the next section. For now, 
let us focus on the geometry related to the trajectory of 
the semihard parton. 

At impact parameter b the initial system is almond 
shaped with width and height being 

w^RA-h/2 , h^[R\~{h/2fY'^ , (1) 

where Ra is the nuclear radius. In the following we 
shall use dimensionless length variables by normalizing 
all lengths by Ra, so ([1|) is to appear as 

w = 1-6/2 , h^{l~b'/^f'^ . (2) 

An ellipse evolving from such a system is to be described 
by the equation 

The initial configuration corresponds to m = 1, so we use 
{xi^yi) to denote the coordinates of that ellipse. Since 
a semihard parton created at P is at a distance t on the 
straight line at angle (ps from the boundary, there is a 
short time interval for the parton to go that distance to 
reach the boundary. We ignore the small expansion that 
the system may undergo during the transit time, since P 
is not far from the boundary. The time when hadroniza- 
tion occurs is much later, when the medium density is 
lower, and u larger. Compared to that time scale, the 
system is almost static, while the parton traverses the 
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medium. This static approximation renders the determi- 
nation of t much easier. 

The point at which the trajectory intersects the bound- 
ary is given by 



Xi — xq + t COS , Ui = Vo + t sin ( 



(4) 



Using Eq. ^ for {xi,yi) at w = 1, we can solve for t, 
getting 



t 



where 



^ = I — cos < 
w 



■ smf 



B 



C = \- {x^/wf - (yo/hf 



(5) 



(6) 



(7) 



(8) 



At any point (x, y) the local flow direction is specified 
by the gradient of u{x, y) even for m < 1. The azimuthal 
angle of that flow direction will be denoted by ip{x,y), 
whose value is 



V'(a;, y) — tan 



'uP'y 
h^x 



(9) 



Since for any b the angle 'ijj{x, y) can vary from — 7r/2 to 
-|-7r/2 depending on the position of the point (x,?/) with 
a; > 0, there can always be a semihard-parton trajectory 
with a (j>s that coincides with tl>{x,y). But it also means 
that there is a broad range of possibilities where (ps dif- 
fers from -0. To emphasize this difference is the main 
characteristic of this model. 

In the following section we shall have dynamical reason 
to follow the flow direction starting from any point (x, y) 
on the parton trajectory along (pg. Let t' denote the 
distance from (x, y) to the surface in the direction 'ip{x, y); 
it can be calculated in the same way as in Eqs. (H])-® 
and the equation for t' is the same as in Eq. ^ except 
for the replacement of 4>s by ip{x,y). 

We shall take the local density D{x, y) in the transverse 
plane to be as described in the Glauter model for AB 
collision 



9AB \ b,s) = Ta{s) 

-Tb(|s-6|) 



(10) 



where T!a(s) is the thickness function normalized to A^ 
i.e. 

Ta[s)^ A j dz p{s,z) , JdhTA{s)=A. (11) 



p is the nuclear density normalized to 1, for which we 
adopt the Woods-Saxon form 



Pir) = Po 



1 + e('-'-o)/? 



(12) 



where ro — 6.45 fm and ^ — 0.55 fm. We shall take the 
effective nuclear radius to be Ra = 7 fm, which is also 
used to scale all length variables, so that when scaled 
the corresponding variables are tq ~ 0.92 and ^ — 0.08. 
The scaled mean density po is then 0.285, and a being 
the inelastic nucleon-nucleon cross section (taken to be 
40 mb) becomes 0.082. For a point {x, y) in the right-half 
almond region we have 



= (x + 6/2)2 + ^ 



(13) 



so the longitudinal lengths of A and B at that point are 
LA.B{x,y) = — dzp{s,z) , (14) 

Po J-ZA,B 

where 



zl = l-s\ 4 = i-|s-5|2 



Relating Eq. ([TT]) to ([HI) we have 
aTAis) ^ ujLAix,y), lo ^ 



4.6 



(15) 



(16) 



where A = 197 has been used. Thus apart from an over- 
all normalization constant, the local density in the trans- 
verse plane is 



D{x,y) = LA{x,y) 
+LB{x,y) 



1 



'U}LB{x,y) 



(17) 



The application of D[x, y) below will not rely on its ab- 
solute magnitude. 



III. RIDGE FORMATION 

Having described the geometry related to the semi- 
hard parton, we proceed now to the consideration of 
hadronization and ridge formation. Since there is no the- 
oretical framework in which one can reliably treat how 
the semihard parton interact with the medium and how 
the energy loss is converted to ridge particles, we pro- 
pose a model that describes the subprocesses in terms 
of relevant distributions with parameters to be deter- 
mined phenomenologically. The conversion of lost energy 
to ridge particles is described by a correlation function, 
which turns out to be central to the phenomenology of 
0s dependence. 

Let a semihard scattering occur at P in the transverse 
plane with a scattered parton moving at angle (pg. The 
distance from P to the boundary along the straight-line 
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trajectory is t. The probability P{xo,yo, t) of detecting a 
parton emerging from the medium is the product of the 
probabihty of producing a semihard parton at {xo,yo), 
which is proportional to the product of the longitudinal 
lengths at that point, La{xo, ?/o)-^s(a^o, ?/o)j and the sur- 
vival probability S(t), i.e., 



P{xQ,yQ,t) (X LAixQ,yQ)LB{xo,yo)S{t) 



(18) 



The proportionality factor that depends on the semihard 
scattering cross section will be canceled when we calcu- 
late the per-trigger yield, so it is not important to have 
it specified here. We assume that S{t) has an exponen- 
tial dependence on t due to the opaqueness of the dense 
medium 



S{t) = exp[-</T(a;o,yo)] 



(19) 



where r (0:0,2/0) should depend on the density along the 
trajectory. To implement the calculation in a manageable 
way, we first evaluate the density function D{xq, j/q) at all 
grid points within the geometric region of interest, and 
then between any given initial point P and the exit point 
we evaluate the local density along the corresponding t- 
segment by means of 2D interpolation among neighboring 
grid points. Thus we write 



T(a;o,2/o) = to/d{xo,yo) , 



(20) 



where is a free parameter and d{xo,yo) is the relative 
density 



dixo,yo)^Dixo,yo)/D{0,0) 



(21) 



with D{0,0) being the density at the center of the over- 
lap X = y = 0. Since only semihard partons created 
near the surface are likely to lead to a trigger particle, 
^0 is expected to be small, so we shall use just one such 
parameter for all trajectories at all centralities. Thus, 
if t is large compared to T{xo,yo), the parton would be 
absorbed by the medium, and cannot leave it to form 
a trigger particle or any structure above the bulk back- 
ground. In the numerical computation we cut off t at 
2T{xQ,yQ). 

To initiate a description of the soft interactions that 
generate the ridge, let us first use to denote the dis- 
tance from P along that trajectory so that ^ = 1 is at the 
boundary. With (xo,?;o) being the coordinates of P, the 
coordinates (a;^,y^) at ^ are 



x^ — Xq + ^tCOS( 



% = yo + £.t sin ( 



(22) 



The probability that the semihard parton emits a soft 
parton at ^ is proportional to I?(x^, y^). We cannot spec- 
ify in more detail the nature of the soft emission in the 
absence of a quantitative description of the soft process. 
We assume at the qualitative level that some gluons are 
emitted by the semihard parton that do not significantly 
alter the straight-line trajectory of the parton. Such soft 
emissions occur at successive points along the path. If 



they augment one another coherently, a significant effect 
may accumulate and lead to observable consequences on 
hadronization. 

Usual perturbative theory applied to the study of en- 
ergy loss of hard partons traversing dense medium is not 
concerned with what happens to the medium. Our con- 
cern here is the opposite. The gluons radiated by a semi- 
hard parton are absorbed by the medium, thereby en- 
hancing the thermal motion of the medium partons in 
the vicinity of the trajectory. Since the medium expands, 
those thermal partons flow collectively and carry the ex- 
tra energy gained along the flow whose direction can be 
determined locally. We are interested in the thermal par- 
tons because of our focus on ridge particles, which have 
the empirical characeristics that 

(a) their p^^°'^ distribution is exponential with slope 

harder than that of the bulk background, 

(b) their yield increases with centrality, and 

(c) the B/M ratio in the ridge is similar to that in the 

inclusive distribution. 

All three combined strongly suggest that the ridge par- 
ticles are formed by recombination of enhanced thermal 
partons jl6|] . The semihard parton that emerges from the 
medium leads to the creation of the trigger particle, but 
plays no direct role in the formation of the ridge. The 
above discussion refers to the average over all triggered 
events. 

We now direct our attention to the dependence on (j)^ . 
In the discussion above two directions are emphasized: 
one is the azimuthal angle (j)s of the trajectory; the other 
is the flow direction, denoted by ij: in the preceding sec- 
tion. The former refers to the semihard parton, while the 
latter refers to the movement of the local medium that 
carries the soft partons in a direction that may or may 
not differ from 0s- If V'l^:??/) is approximately equal to 
(j)s for most of the points {x^y) along the trajectory of 
the semihard parton, then the thermal partons enhanced 
by successive soft emissions are carried by the flow along 
in the same direction; the effects reinforce one another 
and lead to the formation of a ridge in a narrow cone. 
On the other hand, if the two directions are orthogonal, 
then the soft partons emitted from the various points 
along the trajectory are dispersed over a range of sur- 
face area, so their hadronization leads to no pronounced 
effect. These extreme possibilities suggest a correlation 
function between 0s and which we assume to have the 
Gaussian form 



C{x,y,(j)s) = exp 



(0s - i^{x,y)Y 
2A 



(23) 



where the width-squared A is a parameter to be deter- 
mined. This is a phenomenological formula that cannot 
be derived from flrst principles, but has sound physical 
basis and will play a central role in our model. 

To provide a pictorial impression of the correlation be- 
tween 0s and tp that can affect the ridge formation, we 
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(a) 



(b) 



FIG. 1: (Color online) Illustrations of the relationship be- 
tween the trigger directions (j)s in (red) arrows and the flow 
directions t/i in thick (green) arrows for noncentral collision, 
(a) Semihard partons at 0s = originated from 3 different 
points P where only the middle one has matching ifis and ip 
that lead to strong ridge, while in (b) for ~ 70° only the 
upper one has matching angles, leading to stronger ridge than 
in the two lower non-matching cases, but it is weaker than the 
middle one in (a) because of lower local density at the tip of 
the ellipse. 



show in Fig. 1 some examples of possible values of those 
angles. It is important to recognize that the point P of 
creation of the semihard parton can vary over all points 
in the ellipse and that the possible mismatch between 
those angles depends on P. In both panels of Fig. 1 we 
show three illustrative points of P. In (a) we set (ps = 0° 
shown by the thin arrows in red, and in (b) (j)^ ^ 70°. 
The thick arrows indicate the flow directions V' normal 
to the surface. Evidently, the middle P in (a) and the 
upper one in (b) result in matching 0s and ip, while the 
others do not. Ridge formation is stronger for the match- 
ing pairs than the mismatched ones. The ridge in (b) is 
weaker than that in (a) because the density near the top 
of the ellipse is lower than that in the middle. In essence, 
the origin of the 0s dependence of the ridge is depicted 
in Fig. 1. 

For every point (x, y) on the trajectory, the flow di- 
rection ^'{x, y) specifies only the average direction of the 
ridge hadrons, since there are statistical fluctuations, the 
magnitude of which depends on how far {x,y) is away 
from the surface along the direction ^{x,y). That dis- 
tance is t', noted already in the preceding section. We 
introduce another distribution to describe the fluctua- 
tion of the azimuthal angle of a ridge particle from the 
average flow direction 



T{x,y,< 



exp 



{(j>- ip{x,y)f 
27t' 



(24) 



width. Clearly, the farther the emission point is away 
from the surface, the wider fluctuates from ip[x,y). 

We now assemble all the pieces that describe the var- 
ious subprocesses and write the ridge-particle distribu- 
tion i?(0, 0s, 2^0, J/o) due to a semihard parton created at 
{xQ^yo) moving at 0s as a product of all the factors, in- 
tegrated along the trajectory: 

R{(l),(j>s,xo,yo) = NP{xo,yo,t)t 

X [ d£_D{x^,y^)C{x^,y^,(t)s)r{x^,y^,(t)) , (25) 
Jo 

where N is an overall normalization constant. The vari- 
ables i, x^ and y^ all depend implicitly on the initial 
coordinates {xo,yo)- For the observed distribution it is 
necessary to integrate over all (a;o,yo). Not every semi- 
hard parton included in that integration gets out of the 
medium to generate a particle that triggers the event. 
The ridge distribution per trigger is therefore normalized 
by the probability of the ridge-generating parton emerg- 
ing from the medium 



i?(0,0s 



/ dxQdyoR{(f>, 0s,xo,yo) 
J dxodyQP{xo,yo,t) 



(26) 



In experimental analysis of the data positive and negative 
values of 0s are combined to increase the statistics in the 
determination oi dN/dA(p. We shall flnd, however, in our 
calculations interesting details of R{<p, 0s) that depend on 
the sign of 0s, as will be discussed below. 

In Eq. (|26p the integration of (xq, yo) is over the entire 
initial ellipse specified by Eq. ^ with u = 1. For —7t/2 < 
(j)s < 7r/2, the contribution comes mainly from the right- 
half eUipse, but not entirely. At large \yo\ and xq < 
it is possible for a semihard parton to be emitted with 
< |0s| < 7r/2 and survive to trigger an event. Thus we 
shall integrate over both positive and negative values of 
xo and yo. 

Equations (l25l) and (|26p contain the essential elements 
that affect the ridge distribution so long as the depen- 
dence on p^'® and p^^°'^ is not brought into the open. 
Such dependencies enter into the parameters contained 
in the various factors, especially S{t) and C{x^,y^, 4>s)', 
however, since we shall not vary the ranges of p^^'^ and 
passoc jj-^ ^Yiis paper, but just fix them at the narrow 

range of experimental values 3 < p^'^ < 4 GeV/c and 
1.5 < p^^°'^ < 2.0 GeV/c, we leave the pt dependence 
implicit, and proceed to the confrontation with real data. 



DEPENDENCE OF RIDGE YIELD 



where the degree of fluctuation is specifled by jt'; t' is 
the height of the cone of fluctuation and 7 specifles the 



IV. 



Apart from the overall normalization N in Eq. (|25p . 
there are three essential parameters in our model: to in 
S(t), A in C{x,y,4>s), and 7 in T{x^y,4>). They quantify 
three independent characteristics of the ridge formation 
process: survivability, azimuthal correlation, and fluctu- 
ation, respectively. 



6 



Since the shape of the near-side correlation is insen- 
sitive to the integrated yield for |A(/)| < 1 has been 
analyzed for the ridge part and shown as a function of 
(j)s for IA77I < 0.7 and for two centrality bins fl^. Thus 
we define accordingly 



for both central and mid-central collisions, but especially 
for 0-5% centrality. It is nontrivial that the decrease with 
increasing 0s can be so different for the two centrality 
cases and then they both flatten out above 60°. 



R{4>, (t>s) 



(27) 



as the ridge yield per trigger for any impact parameter h. 
Carrying out the above integration over i?(0, 0s, xq, ?/o)j 
where according to Eq. ([^S]) only the factor T{x^,y^, cf)) 
depends on 0, we obtain 



del) R{4>, 0s, xo, yo) = NP{xo, yo, t)t 



X d^ D{x^,y^)C{x^,y^,(t)s)G{x^,y^,(j)s), (28) 
/o 



where 



Gix^,y^, 



,+1 



exp 



(0-0(a:g,y^))^ 
27t' 



(29) 

The last integral is not sensitive to ^, 0s or 7, so 
G{x(^,y^, acts as a modifier of the overall normaliza- 
tion N. The parameter 7 will be determined in the next 
section by the differential correlation that depends on 0, 
but here we focus first on the integrated yield, using 7=1 
that will be shown to be the final value. The point is that 
the yield i^(0s) depends mainly on two parameters: to 
and A. 

The yield data are shown in Fig. 2 for (a) top 5% and 
(b) 20-60% with 0s in 6 segments ranging from to 90°. 
The mild dependence on 0s in (a) is not surprising for 
central collisions, but the precipitous decrease in (b) is 
striking. The three main features of Fig. 2, namely, the 
relative normalization between the yields of the two cen- 
tralities, the different rates of decrease with 0s, and the 
flattening out at large 0s , are correlated, and we fit them 
by varying the two parameters to and A. The yield for 
the top 5% is calculated by taking the average of 6 = 0.3 
and 0.4. The yield for the mid-central case is obtained by 
averaging over b — 1 and 1.3, corresponding to the wider 
experimental range of centrality 20-60%. The solid lines 
in Fig. 2 represent the best fit we can achieve with values 



to 



0.2, 



A = 0.11 



(30) 



In Fig. 2(a) the height of the solid line is adjusted to 
fit the data point at the lowest 0s by varying the nor- 
malization factor N in Eq. (pS)) . TV encapsules all the 
uncalculable effects of the soft processes involved in the 
ridge formation, and is not essential to the study of the 
0s dependence. However, once it is fixed at A'' = 0.085 by 
the top 5% data, the normalization for the mid-central 
collisions in Fig. 2(b) is no longer adjustable. Our re- 
sults reproduce the characteristics of the data very well 
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FIG. 2: (Color online) Dependence of ridge yield on 0s for (a) 
top 5% and (b) 20-60%. Data are from Ref. [Hi. The solid 
lines are the results of calculation in CEM. 



The values of the parameters in Eq. (|30p are very rea- 
sonable. First, to = 0.2 implies that only semihard par- 
tons created on the average at Ra/4: from the surface or 
closer get out to form trigger and ridge particles. Second, 
A = 0.11 implies that the width Cc = VX of the corre- 
lation function is 0.34 rad, a value that is significantly 
smaller, as it should, than A0 ~ 1, where the wings of 
the ridge vanish. The most important piece of physics 
that the phenomenology has revealed to us is contained 
in the correlation function C(a;,y, 0s) given in Eq. (|23p . 
It implies that a strong ridge can be formed only if the 
local flow direction ip{x,y) is within a cone of about 20° 
from the direction of the semihard parton 0s . The rapid 
decrease of the yield with 0s seen in Fig. 2(b) is a direct 
consequence of the mismatch between the two directions. 
The degree of that mismatch depends on the initial point 
(xq, yo) for any fixed 0s. As depicted in Fig. 1, if 0s ~ 0, 
the mismatch is larger at large |yo| than at small |yo|j 
whereas if 0s ~ 7r/2, the opposite is true. In the former 
case matching condition occurs (i.e., when 0s ~ i/i ~ 0) 
at small |yo|, where the density is high, thus enhancing 
soft emission; in the latter case (when 0s ~ "0 w 7r/2) 
it occurs at large yo, where the density is low, thus sup- 
pressing ridge formation. That is why the ridge yield 
decreases with increasing 0s, and does so more rapidly 
in noncentral collisions because the density varies more 
significantly near the surface of the overlap in that case. 
The leveling-off of the descent at high 0s in both cen- 
trality cases is due to the contribution from the partons 
created in the left-half ellipse at high yo- 
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V. Acj> DISTRIBUTION OF THE RIDGE 

Having reproduced the data on the (f>s dependence of 
the ridge yield, we now consider the Acj) dependence of 
the ridge particles. To that end we study the behavior of 
Ricj), 0s) in Eq- (|26|) without integrating over as done in 
Eq. (gt]). The data in Ref. show that the near-side 
peaks in {l/Ntrig)dN/dA(j>, where A(j> = (j> — (j)s, have 
rather similar shape for different bins of <j)s and p^^""^. 
For definiteness, we focus on the bin for 15° < <j>s < 30° 
at 20-60% centrality, as shown in Fig. 3(a). They in- 
clude both jet and ridge components. The dashed line is 
the calculated result for R{(j),(f)s) at </)s = 22° (left curve) 
and the dashed-dotted line is for ips = —22° (right curve), 
plotted as functions oi <j) — 4>s- Both are calculated for 
7=1. Due to the symmetry in the problem the two 
curves are mirror reflections of each other about Acf) = 0. 
Since the data on the A(f) distribution include the con- 
tributions from both positive and negative values of (j)s, 
we average the two curves and obtain the solid line for 
|(/)s| = 22°. The difference between the data and the cal- 
culated R{A(j), \(j)s\) is the jet contribution, J {Acj), \(j)s\)- 
In the visual presentation of Ref. [l^, not included in 
the proceedings, the data for the ridge distribution are 
shown; the segment for 15° < (j)s < 30° at 20-60% cen- 
trality is reproduced in Fig. 3(b). The solid line in that 
plot is our calculated result for i?(0, (j)s) with 7=1 [iden- 
tical to that in Fig. 3(a)], showing good agreement with 
the data. 
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FIG. 3: (Color online) The data are A(j) distributions from 
[ll] for 15 < 0s < 30° at 20-60% centrality for (a) the sum of 
jet and ridge and (b) ridge alone. The curves are all calculated 
in the CEM for the ridge distributions only with 7 = 1. The 
dashed and dashed-dotted lines are left- and right-shifted for 
(j>s = ±22°, respectively. The solid lines are the average over 
the two signs of (j)s- 

To see the sensitivity of our result to the value of 7, we 
show in Fig. 4 two curves for i?(A0, \(j)s\) for 7 = 1 (solid) 
and 7 = 2 (dashed). In view of the large error bars both 
are acceptable, although the solid line agrees better with 
the data on both sides of A(j) = 0. Hereafter, we shall 



regard 7 = 1 as the value determined from fitting the 
correlation data. 



0.12 - 




FIG. 4: (Color online) Same as Fig. 2(b) with the solid line 
for 7 = 1 and dashed line for 7 = 2. 



The asymmetric result for positive and negative 4>s 
shown in Fig. 3(a) is significant. For positive (ps, the shift 
of the peak of i?(A0,0s = 22°) to the left of Acj) ^ 
shown in Fig. 3(a) can be studied for other values of 
(j)s also. We show in Fig. 5 our calculated results for 
R{(l),(t>s) at <j>s = 7°, 22°, 37°, 52°, 67° and 82°. The left- 
ward displacements from A(f) = reach their maximum 
at (f>s w 37° where the magnitude of the shift is approxi- 
mately 10°. We regard this asymmetry of the ridge dis- 
tributions as an indication of some important aspect of 
the physics, which we shall in the next section quantify 
in a way that may be easier for experiments to verify. 
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FIG. 5: (Color online) The ridge distributions for various 
positive values of (j>s- 
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VI. INSIDE-OUTSIDE ASYMMETRY 
FUNCTION 

The shifting of the maxima of the ridge distributions 
seen in Figs. 3(a) and 5 reveals some property of the 
correlation between semihard parton and the ridge par- 
ticles, although the phenomenon is not seen directly in 
the correlation data that collect both positive and nega- 
tive values of (t)s. The ridge yields in Fig. 2 reveal only 
the gross effect of the correlation because of the integra- 
tion over (p in Eq. (j27p . but not some detailed proper- 
ties. To have a qualitative understanding of the origin 
of the shift, let us consider the case (j)s = 7r/4, for in- 
stance, bearing in mind the difference between (jig and 
ip{x,y). For most values of {x,y) along the parton tra- 
jectories, the flow directions normal to the surface have 
ij}{x,y) < 7r/4. However, for (xQ^yo) in a small region 
near the top of the ellipse, ip{x,y) is > 7r/4, but that 
is where the density is low. Thus in most of the regions 
where the semihard partons are produced with 0s — 7r/4, 
the ridge hadrons have (j) directions mostly smaller than 
the parton direction. We therefore expect an asymmetry 
between (p < (pg and 4> > 4>s. To quantify these proper- 
ties in a way accessible to experimental analysis, we now 
introduce a measure that we refer to as the inside-outside 
asymmetry. 

Restricting (ps to the range < (ps < 7r/2 we define 



Y+{ct>s) 



d(j)R{(j), (ps), 



,+1 



d^R{(t>,(t>s) 



(31a) 
(31b) 



But for — 7r/2 < 0s < wc define them in the opposite 
way 



Y( 



d(pR{(p, 0s), 



d(f>R{^,^s) 



(32a) 



(32b) 



Thus y+ may be thought of as being mostly inside (i.e., 
on the in-plane side of 0t), while Y- may be thought 
of as being mostly outside (i.e., on the out-of-plane side 
of 0t)- Clearly, when 0s = 0, it is necessary that 1"+ 
equals Y_ to be consistent, and of course they are equal 
by reflection symmetry across the horizontal axis. Also, 
at 0s = ±7r/2 reflection symmetry across the vertical axis 
requires 1+ = Y^. We define the asymmetry function to 
be 



y+(0s)-y-(0s) 
r+(0s) + r_(0s) 



(33) 



where the denominator is just Y{(ps) as defined in Eq. 
(I?7)l . It follows from reflection symmetries that A{0) = 
and A{±tt/2) = 0. How A(0s) varies between the two 



extremes at 0s = and tt/2 depends on b. Experimen- 
tally, the two definitions of y±(0s) in Eqs. (27) and (28) 
are necessary in order that the data in the two quadrants 
of 0s may be combined to increase statistics, but theo- 
retically, one is the reflection of the other, so the study 
of the sector < 0s < 7r/2 is sufficient. 
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FIG. 6: (Color online) The asymmetry function y4(0s 
0-5% (dashed) and 20-60% (solid). 



for 



In Fig. 6 we show the asymmetry function for the two 
centrality bins. It is evident that for the entire range 
of 0s the asymmetry A(0s) for mid-central collisions is 
more than twice larger than that for central collisions. 
Of course, a.t b — there should be no asymmetry. The 
rapid growth of A(0s) with increasing 6 is a striking fea- 
ture of the effects of the mismatch between 0s and V' 
that complements the 0s dependence of the ridge yield. 
Experimental verification of our prediction of the proper- 
ties of A(0s) would lend further support to the correlated 
emission model. 



VII. CENTRALITY DEPENDENCE 

We have investigated more thoroughly the centrality 
dependence of the yield per trigger, and show in Fig. 7(a) 
our result for Y(b, <ps) versus impact parameter b for vari- 
ous values of <ps- The monotonic decrease of Y{b, <ps) with 
increasing b occurs only for 0s > 40°. At lower 0s it in- 
creases with b initially before decreasing. Such a peaking 
is an unexpected finding. However, after averaging over 
all 0s the bump in the per-trigger yield disappears. For 
the averaging process we first note that the denominator 
of i?(0, 0s) defined in Eq. (|26p has implicit dependence 
on 0s because t given in Eq. ([5]) does. Thus the average 
yield is obtained by averaging both the numerator and 
denominator of Yib, (ps) over 0s separately, i.e.. 



Y(b) = 



/c 



-it/2 



jr^i dc/) J dxodyoR{(l), 0s, 2:0, yo) 



Jq^^ d(t>s J dxodyoP{xQ,yo, t) 



(34) 
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where the dependencies on b in the integrands are im- 
pUcit. It is this Y{b) that is plotted in Fig. 7(b). When 
plotted against A^part the average yield increases mono- 
tonically, as shown in Fig. 7(c). Since the magnitude 
of the per-trigger yield depends sensitively on the cut in 
^assoc^ we do not show in that figure any data having cuts 
different from 1.5 < p^^°'^ < 2.0 GeV/c. The two points 
in that figure are determined from the data in Ref. p^ . 
using the number of triggers in each (f>s bin provided by 
the group that performed the analysis. 
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FIG. 7: (Color online) (a) Ridge yield per trigger vs impact 
parameter for 5 values of (ps, (b) Y{b), averaged over all (j>s, 
vs impact parameter, and (c) average yield vs A^'part ■ The two 
points in (c) are determined from the data in Fig. 2(a) and 
(b). 

It is of interest to remark on the significance of the 
peak of Y{b,(j)s) at mid-6 for small (jig in Fig. 7(a). The 
peak suggests that at small 0^ the yield can be higher 
in noncentral collisions than in central collisions. That 
is a theoretical result that has phenomenological support 
as can be seen in Fig. 2, where the data point at the 
lowest (f>s at 7° in panel (b) is higher than that in panel 
(a). Our calculated result in Fig. 2(b) for Y{(j)s), after 
integrating over the range of centrality 20-60%, shows 
a slower descent with (pg compared to the data, and at 
(j>s — 7° it is ~ 10% lower than the data point. It means 
that the peaking of our theoretical result on Y{b,4>s) in 
Fig. 7(a), though remarkable in its existence, is not high 
enough. The most likely source of the discrepancy lies 
in the value of A for mid-values of 5, since a decrease 
of A for 1 < b < 1.3 can produce a steeper descend in 
(j)s. However, we prefer at this stage of our investigation 
not to introduce extra freedom to achieve a better fit 
of the data. In our view the dynamical origin of the 
observed phenomenon at low (jjs is more important than 
the numerical accuracy of our result. As we now proceed 
to describe, our model provides an explanation of the 
small (jjs behavior. 

The geometrical difference between mid-central and 



central collisions is, of course, the extra overlap space 
in the latter case. In that extra space the deep inte- 
rior docs not make any contribution to the ridge, since 
the semihard partons originated there cannot survive on 
their way to the surface. The only regions where the 
partons moving at small (ps can get out are at the top 
and bottom of the ellipse; i.e., those partons form tan- 
gential jets. However, the flow vectors in those regions 
are essentially normal to the parton trajectories. Thus 
the correlation function C(x, y, 4>s) suppresses the forma- 
tion of ridge particles for partons originating in the extra 
space. Yet the number of trigger particles is not sup- 
pressed because of the tangential jets. Consequently, the 
yield/trigger decreases with decreasing b. Note that the 
argument does not work if (pa is large because then 'ip{x,y) 
would be close to cps and there would be no suppression 
due to C{x, y, 4>s)- 

Support for our assertion that ridge formation by tan- 
gential jets is suppressed can be found in 3-particle cor- 
relation in dijet events where the two trigger jets have 
nearly equal and opposite momenta. It has been reported 
that preliminary results from STAR indicate the lack of 
ridge structure in events with one trigger at > 5 
GeV/c, another at pt > 4 GcV/c and associated par- 
ticles at > 1-5 GeV/c in AuAu collisions at any cen- 
trality [2l| . From the iVp^^^ dependence of the number of 
dijets triggered, it is suggested that those events are due 
to tangential jets generated near the surface. The anal- 
ysis in Ref. !2l!| makes no reference to (j)s and integrates 
over all (pT- Clearly, tangential jets can be directed at 
any azimuthal angle for any A^part- Our explanation of 
the phenomenon is that for tangential jets the flow direc- 
tion is normal to the semihard parton direction at any <ps . 
We predict that as the momentum difference of the two 
trigger jets increases, corresponding to their becoming 
less tangential, the mismatch between (ps and ip becomes 
less severe, and a ridge should begin to develop. 

A phenomenon related to the dip at small b has been 
observed previously by PHENIX in the dependence of 
the yield on jVpart for baryon triggered jets (without ridge 
separation) [22| . The data show a dip of the yield at high 
A'^part- Such behavior of centrality dependence has never 
been explained before. Our study here cannot address 
the issue of species dependence of the trigger, and the 
data of Ref. [2^ cannot be used to address the issue of 
ridge yield. Nevertheless, our result in Fig. 7(a) holds 
the intriguing possibility that our description of ridge 
formation may contain the germ of an explanation of the 
dip phenomenon. 



VIII. CONCLUSION 

Much theoretical attention has been given to the mod- 
ification of the jet behavior at very high pt due to the 
propagation of a hard parton through dense medium. 
The focus has been on what happens to that hard parton. 
Our interest in this paper has been on what happens to 
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the medium when a semihard parton passes through it. 
It is the complement to the problem at high pt. Both 
are important, since jet-medium interaction acts in both 
ways. The effect of semihard partons on the medium 
is difficult to calculate, so a phenomcnological model is 
needed to relate different existing empirical facts and to 
generate a coherent picture of what the important sub- 
processes are. 

Ridge formation provides crucial evidence that the 
medium's response is correlated to the jet direction. Be- 
cause of the time lag in hadron production, local hydro- 
dynamical flow of the medium can influence the direction 
in which the lost energy goes and where the ridge is to be 
formed. We have attempted to capture these aspects of 
the dynamics by use of some phenomenological functions 
in the CEM. What we have learned is that ridge forma- 
tion can be strong only if the direction of the flow that 
carries the enhanced thermal partons is within a narrow 
cone of about 20° around the semihard-parton direction. 
Within that cone the effects of successive soft emissions 
by the parton reinforce one another to develop a ridge. 
When a parton's trajectory is normal to the flow direc- 
tion, the effect of energy loss is spread out over a wide 
spatial domain and the detection of any ridge formation 
is suppressed. 

We have shown that the CEM reproduces the <j)s de- 
pendence of the ridge yield. Signiflcant insight has been 
gained by examining the details of the ridge character- 
istics. As a result, we have made predictions on certain 
properties of the ridges that can be measured. One is 
on the inside-outside asymmetry behavior relative to the 
trigger direction. Another is on the centrality depen- 
dence of the ridge yield. Still another is on the lack of 
ridge structure in dijet events. Verification of our predic- 
tions will lend additional support to our interpretation 
of the medium response to jets created in heavy-ion col- 
lisions. 

Although hadronization by recombination has not been 
used explicitly in our calculation, it is a subprocess that 
is implied in our modeling. The conversion of the en- 
ergy lost by the semihard parton to the enhancement of 
thermal partons cannot be described quantitatively, so 
the angular distribution of the hadrons formed by the 



recombination of the thermal partons cannot be deter- 
mined in the absence of calculable distribution of the soft 
partons. The dependence of the ridge has been stud- 
ied phcnomenologically in Ref. Q , when the dependence 
on 0s is not an issue. Here we have gone beyond that 
first attempt and have focused on the (f>s dependence. 
Our use of the correlation function C{x,y,4's) is based 
on the identification of the average hadronic direction 
with the local flow direction tp{x,y) of the soft partons 
by virtue of thermal-thermal recombination. One may 
argue that the identification of the average hadronic di- 
rection with the fiow direction is also an attribute of other 
schemes of hadronization, such as Cooper- Frye and local 
parton-hadron duality. However, it is only in thermal- 
thermal recombination that the B/M ratio can be as 
high as observed p3 , [l^ . If ridge formation were domi- 
nated by thermal-shower recombination, we would have 
had to calculate the shower parton distribution, as done 
in Ref. [23.] for V2 at intermediate px, and the hadron 
direction in the ridge would not be simply determined 
by the flow direction of the thermal partons. Clearly, if 
ridges were formed by means of fragmentation of hard 
or semihard partons, there would be no suppression due 
to a mismatch of parton and flow directions. As a con- 
sequence, there would not be any dependence of the 
ridge formation as observed in the experiment. 

We have not investigated in this paper the dependence 
on Pt for either the trigger or the ridge particles. Neither 
have we considered the hadron species of those particles. 
For all Pt < 6 GeV/c recombination is the only viable 
mechanism for hadronization [16]. To relate the various 
dependencies (i.e., on py, 0, 0s, 77, A^part, hadron species 
in trigger, ridge, and jet) is still a larger problem that 
remains to be studied. 
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